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Abstract. We expand the Tsallis distribution in a Taylor series of powers of (g — 1), where q is the Tsallis 
parameter, assuming q is very close to 1. This helps in studying the degree of deviation of transverse 
momentum spectra and other thermodynamic quantities from a thermalized Boltzmann distribution. Af¬ 
ter checking thermodynamic consistency, we provide analytical results for the Tsallis distribution in the 
presence of collective flow up to the first order of (g — 1). The formulae are compared with the experimental 
data. 
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1 Introduction 

It is now a standard practice to use the Tsallis distri¬ 
bution [T] for describing the transverse momentum dis¬ 
tributions at high energies. This has been pioneered by 
the PHENIX and STAR collaborations [11[31I3] at the Rel¬ 
ativistic Heavy Ion Collider (RHIC) at BNL and by the 
ALICE, ATLAS and CMS collaborations [5l[6l[7ll8l[^[T0] at 
the Large Hadron Collider (LHC) at CERN. The Tsallis 
distribution is successful in explaining the experimental 
transverse momentum distribution, longitudinal momen¬ 
tum fraction distribution as well as rapidity distribution 
of hadrons off the e+e“ as well as p — p collisions 
[niiiiiiioreiin]. The form of the Tsallis distribution used 
in this paper has been described in detail previously [151 
muiniEi] and has the advantage of being thermodynami¬ 
cally consistent. There is clear evidence for a mild energy 
dependence of the parameters q and T m- Also, initially 
there were indications that the values obtained for the pa¬ 
rameters q and T were consistent with each other for dif¬ 
ferent particle species dZldl]. Different conclusions have 
been reached in the literature [22l[28ll24ll2R[ . albeit using 
slightly different formalisms and approaches, and a more 
detailed analysis is still outstanding to prove this beyond 
doubt. 


The relevant thermodynamic quantities can be written 
as integrals over the following distribution function: 

1 

q-1 

( 1 ) 


/ = 


l + (9-l) 


E — p 


It can be shown m that the entropy, S, particle number, 
iV, energy density, e, and the pressure, P. 
are given by, 


(tp /I. 

( 2 ) 
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where V is the volume and g is the degeneracy factor. 

The function appearing in Eq. ([^ is often referred to 
as q-logarithm and is defined by 


2 Review of the Main Ingredients of the 
Model 

For completeness we recall here the main ingredients. 
Corresponding Author E-mail: Raghunath.Sahoo@cern.ch 


The first and second laws of thermodynamics lead to 
the following two differential relations: 

de = T ds + p dn, (6) 

dP = s dT + n dp. (7) 

where, s = S/V and n = N/V are the entropy and particle 
number densities, respectively. 
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It is seen that if we use in stead of / to define the 
thermodynamic variables, the above equations satisfy the 
thermodynamic consistency conditions which require that 
the following relations to be satisfied: 


de 

dP 

dP 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 


Eq. ([^, in particular, shows that the variable T ap¬ 
pearing in Eq. Q can indeed be identified as a thermody¬ 
namic temperature and is more than just another param¬ 
eter. It is straightforward to show that these relations are 
indeed satisfied [19]. 

Based on the above expressions the particle distribu¬ 
tion can be rewritten, using variables appropriate for high- 
energy physics as 


dN 

dprdy 


gv 

(27r)2 


PTiriTCOshy 
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mrcoshy — p 
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( 12 ) 


It can be shown that at central rapidity, y = 0, one can 
obtain the transverse momentum distribution in terms of 
the central rapidity density, dN/dy\y^o, as the volume de¬ 
pendence can be replaced by a dependence on dN/dy\y=o 
using 
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gV 
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m — II 
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This leads to the following expression and generalizes 
the expression given in |26j to non-zero values of the chem¬ 
ical potential p (see Appendix for an outline of the 
derivation of Eq. 0 ): 
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where the abbreviations d = m — p and S = 3 — 2q have 
been used. 

In all fits to transverse momentum spectra, the pa¬ 
rameter q turns out to be very close to 1 Ennui- In fact, 
the value of the non-extensive parameter q for high en¬ 
ergy collisions is found to be 1 < q < 1.2 EHEil. In the 
limit where q is exactly 1, Eq. (12) reduces to the standard 


exponential function appearing in the Boltzmann distribu¬ 
tion. It is therefore useful to expand the above expressions 
in a Taylor series in {q — 1) and see how the deviations 
from a Boltzmann distribution develop. Such an expansion 
has been considered previously in EHlUSj- The present pa¬ 
per develops a more systematic analysis than the previous 
ones and considers a slightly different form of the Tsallis 
distribution, having an extra power of q, because it is con¬ 
sistent with basic thermodynamic relations. 

The aim of this paper is to develop a Taylor expansion 
of Eq. (121 in {q — 1) based on the fact that (g — 1) << 1 
(see for example EOl). The conditions of validity of such 
an expansion for pure Tsallis distribution (Eq. Q) is |1 — 
q\E/T < 1. Apart from this, up to first order in (g— l) an 
additional condition |1 — q\{E/TY < 2 must be satisfied 
[28] . The condition of validity for expansion up to order 
(g — 1)^ term will be |1 — gp(A/r)^ < 3. The expansion 
to higher orders has also been considered in |31j in the 
framework of an analysis of quasi-additivity for the Tsallis 
entropy for different subsystems. 

The Taylor expansion is useful as a mathematical tool 
because it breaks the Tsallis distribution in a series of 
(g — 1) containing powers of energy E. Now, the advan¬ 
tage we get is, it will be easier to consistently include the 
effect of flow on the Tsallis distribution just by making a 
substitution E p^Uy,, for a collective four-velocity of 
particles with representative four-momentum p^ Eiiisa- 

There have been earlier attempts to include the effect 
of collectivity in the dynamics of the particles following 
the Tsallis distribution in the form of Tsallis-Blast Wave 
(TBW) formalism Eljj or in Refs. [5^134] . In all the cases, 
an ansatz of fluid four-velocity is taken and energy is re¬ 
placed by the scalar product p^Uy. The inclusion of flow 
inside non-extensive statistics reduces the value of g jSS] 
since some degree of non-extensivity is shared by the dy¬ 
namics. Also, whenever we have an inhomogeneous ther¬ 
modynamic system, with regions having different temper¬ 
atures and exchanging heat with the bigger system, we can 
define an effective temperature which is affected by 
energy transfer only when g 7 ^ I [M1I371I38] . The variation 
of effective temperature with g is seen in |39] . Another im¬ 
portant observed phenomenon like rriT scaling is affected 
by flow. In p — p ov p — A collisions, the particle spectra 





































T. Bhattacharyya et al: Non-Extensive Thermodynamics 


3 


for different hadrons are having the same slope parameter 
T and this phenomenon is known as rriT scaling |40| . Due 
to this scaling behaviour, the particle species cannot be 
identified just by looking at the spectrum. This behaviour 
will be manifested in particle distributions following Eq. 
(12). But, because of the inclusion of collectivity the par- 
tiHe species start having different slope parameters. The 
average collective velocity as well as the mass of parti¬ 
cles will contribute to the slope parameter and hence the 
scaling is broken m- It is possible that q represents the 
joint action of many dynamical factors |42) and it cannot 
be excluded at present that when these are accounted for, 
the factor q — 1 could well become zero. The values of 
the parameters are clearly sensitive to the details of the 
hadronization mechanism gSEl!. 

The TBW formalism considers an implicit dependence 
of fluid rapidity on the velocity; but in the present case 
we consider a cylindrical geometry and velocity comes ex¬ 
plicitly in the calculations. In contrast to the numerical 
treatment in [22j , we provide an analytical formula for the 
Tsallis distribution expanded up to 0{q— 1) in presence of 
flow. In the next section, we derive the Taylor expansion 
for the Tsallis distribution in the series of {q — 1) and we 
compare the order by order deviation from a thermalized 
Boltzmann distribution of transverse momentum spectra 
in hadronic and nuclear collisions. In Sec. 4 we find out 
the expressions for number density (n), pressure (P) and 
energy density (e) for a system with non-zero chemical 
potential (/r). We verify the thermodynamic relationship 
n = dP/d^. In Sec. 5 we derive an analytical formula in¬ 
cluding flow up to ©(q — 1) of a Taylor series expansion 
of the Tsallis distribution and show its application to ex¬ 
perimental data. Finally we summarize the paper in Sec. 
6 . 



3 Momentum Distributions to First and 
Second Order in {q — 1 ) 

Assuming the parameter q to be close to 1, as is the value 
for many cases in high energy physics, the modified Tsal¬ 
lis distribution, /'^, appearing in the expressions for the 
thermodynamic quantities can be expanded in a Taylor 
series in g — 1 with the following result (see appendix 
for a detailed derivation): 


Fig. 1. Fits to the normalized differential yields of as mea¬ 
sured by the ALICE collaboration in p — p collisions at ^/s = 
0.9 TeV [S] fitted with the Tsallis (solid line) and Boltzmann 
distributions (dashed line). Also shown are fits with the Tsallis 
distribution keeping terms to first (dash-dotted line) and sec¬ 
ond order in {q — 1) (dotted line). The lower part of the figure 
shows the difference between model (M) and experiment (E) 
normalized to the model (M) values. 
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dN 

d^p 


= CE 


i + (g-i) 


E — pL 
T 


g 

q-l 


(16) 


where C = gV/{2'KY. Let us use the following notations: 
(g-l)= X] {E-pl)/T = <P axiA l + {q-l)^ = l+x<P = 
f(x). Hence, the expansion of the Tsallis distribution up 
to 0{x’^) can be written as: 








This result can be used for the invariant yield of par¬ 
ticles if it is written in terms of the Tsallis distribution. 


(17) 
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Fig. 2. Fits to the normalized differential yields of tt"*", includ¬ 
ing only results at small values of the transverse momentum, 
as measured by the ALICE collaboration in p — p collisions at 
= 0.9 TeV fitted with the Tsallis (solid line) and Boltz¬ 
mann distributions (dashed line). Also shown are fits with the 
Tsallis distribution keeping terms to first (dash-dotted line) 
and second order in {q — 1) (dotted line). 


Hence one obtains 
dN 


, X 


2 ^2 

-b CE^ — (24 - 20^ -b (18) 



Since, we use a modified form of Tsallis distribution, 
a comparison with similar work will be worthwhile at this 
point. With this aim, we integrate over the rapidity vari¬ 
able to compare the transverse mass spectrum obtained 
from the present approach with that obtained in Ref. |H] : 
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mTdrriT 
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(19) 


Fig. 3. Fits to the normalized differential tt yields as mea¬ 
sured by the ALICE collaboration in (0 — 5)% Pb-Pb collisions 
at .^snn = 2.76 TeV [101 fitted with the Tsallis (solid line) and 
Boltzmann distributions (dashed line). Also shown are fits with 
the Tsallis distribution keeping terms to first (dash-dotted line) 
and second order in {q — 1) (dotted line). The lower part of the 
figure shows the difference between model (M) and experiment 
(E) normalized to the model (M) values. 


fast with increasing y. And so, according to the standard 
practice, we can effectively replace the rapidity integration 
from 0 to ymax to 0 to oo so that the integration yields 
Bessel functions. The first term in the above expression is 
the well-known formula for a thermal source with a Boltz¬ 
mann distribution: 


dN 

PrdpT 


EL 

27r2 


Ki 



( 20 ) 


where iC„s are the modified Bessel’s functions of second 
kind (see Appendix^. In principle, while examining the 
transverse spectra, the rapidity integration should be within 
a maximum value ymax, say. But, owing to the presence 
of the term g-"iTcoshi//T^ integrand drops down very 


In the limit /r = 0, the transverse mass spectrum obtained 
from Ref. [29] (up to 0{q — 1)) does not contain the term 
involving Kq and K 2 . 

The Boltzmann distribution, the pure Tsallis distri¬ 
bution (Eq. (HI) and the expansion of the Tsallis distri- 
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but ion up to the first and second order of {q — 1) (Eq. 
(18)) were used to fit the transverse momentum distribu¬ 
tions obtained by the ALICE collaboration. The results 
are shown in Fig.[^for p-p collisions at ^/s = 0.9 TeV and 
in Fig. I^for Pb-Pb collisions at ^snn = 2.76 TeV. It is 
well known that the Tsallis fits give excellent results for 
p-p collisions but are not very good for Pb-Pb collisions 
as can be seen clearly in Fig. (2) where a large deviation 
at small px is seen. 

It can be seen from the fits to the p—p distribution that 
the successive terms in (q — 1) improve the fits but not in 
a convincing manner. Clearly, it is the best to use the full 
Tsallis distribution and not the series expansion. It might 
turn out of course that the series expansion could be of 
use in a different situation where the comparison with a 
Boltzmann distribution is more relevant. 


The fits in the figures were done using the MINUIT 
package with the following numerical results. 

In Fig. we show fits to the transverse momentum 
distribution of 7 r+ in p-p collisions at 900 GeV. For the 
plain Tsallis distribution (solid line) the parameters were 
obtained as being T = 70.8 MeV, q = 1.1474. The volume 
V was determined as corresponding to a spherical radius 
of 4.81 fm. For the Boltzmann distribution (dashed line) 
the parameters were determined as being T = 150.2 MeV, 
while the radius used to determine the volume was fixed 
at a value of 2.65 fm. For the fit using the Boltzmann dis¬ 
tribution and the first order term in (g — 1 ) (dashed-dotted 
line) the values are T = 138.4 MeV, q = 1.035 while the 
radius is given by 2.80 fm. In the last case corresponding 
to Boltzmann plus first and second orders in (q—l) (dot¬ 
ted line) one has T = 121.2 MeV, q = 1.065 and a radius 
of 3.09 fm. As is well-known and evident, the fit using the 
Tsallis distribution is very good. 

In Fig. we show fits to the small transverse mo¬ 
mentum region of 7 r+ in p-p collisions at 900 GeV. In this 
case the fits using the Boltzmann distribution are fairly 
good. The deviation with the Tsallis distribution becomes 
very prominent only for larger transverse momenta. For 
the plain Tsallis distribution (solid line) the parameters 
were obtained as being T = 70.8 MeV, q = 1.145, simi¬ 
lar to the full range of px- The volume V was determined 
as corresponding to a spherical radius of 4.82 fm. For the 
Boltzmann distribution (dashed line) the parameters were 
determined as being T = 104.9 MeV, i.e. much lower than 
the full range in px while the radius used to determine the 
volume was fixed at a value of 3.61 fm. For the fit using the 
Boltzmann distribution and the first order term in {q— 1) 
(dashed-dotted line) the values are T = 89.9 MeV, q = 
1.07 while the radius is given by 4.03 fm. In the last case 
corresponding to Boltzmann plus first and second orders 
in {q— 1) (dotted line) one has T = 77 MeV, q = 1.11 and 
a radius of 4.54 fm. 

In Fig. [^we show fits to the normalized differential 7 r“ 
yields in (0 — 5)% Pb-Pb collisions at ^snn = 2.76 TeV as 
measured by the ALICE collaboration [10] with the Tsal¬ 
lis (solid line) and Boltzmann distributions (dashed line). 
Also shown are fits with the Tsallis distribution keeping 
terms to first order (dash-dotted line) and second order in 


(<7 — 1) (dotted line). The lower part of the figure shows 
the difference between the Tsallis distribution (M) and ex¬ 
periment (E). It is clear that the best fit is achieved with 
the full Tsallis distribution, whereas, using the Boltzmann 
distribution the description is not good. Successive correc¬ 
tions in (g — 1) improve the description. There is a clear 
deviation at very low transverse momentum (below 0.5 
GeV) and also at higher values above 2.75 GeV. 


4 Thermodynamic Quantities to First Order 

in (g - 1) 

The particle density in Tsallis thermodynamics is given to 
first order in (g — 1) by the following expression: 

n^ + (g-l)n^ (21) 

where is the standard Boltzmann result for the particle 
density: 

= £^e^T^a^K2{a), ( 22 ) 

with a = m/T, and the first order expression in g — 1 is 
given by 

y \—6a^K2{a) — 2a^Ki{a) 

47r^ 

—4a^bK2{a) + 3a^ 7 ^ 3 ( 0 ) -I- a'^K2{a) + a?}P'K2{a) 
-2a?hKi{a)\ . (23) 

In Fig. [^ we show the ratio of the particle density to 

first order in (g — 1) to the full particle density as given 
by the Tsallis distribution, {n^ + (g— l)n^)/n for several 
values of g indicated in the figure as a function of the 
temperature T. 

It can be seen that the expansion in (g — 1) is excellent if 
(g — 1) = 0.01 but rapidly deviates from the full Tsallis 
distribution for larger values of g. Already for (g—1) « 0.1 
the deviations are of the order of 10 % as can be seen from 
Fig. 13 

For comparison we show in Fig. the first order ex¬ 
pansion compared to the Boltzmann expression, (n® -|- 
(g — l)n^)/n^, again as a function of the temperature 
T for several values of the parameter g. In this case the 
deviations are most pronounced for small values of the 
temperature. 

Gorrespondingly, the energy density is obtained as: 

+ (g - l)e^ (24) 


= ^^^2 {ia^K2{a) + a^Ki{a)) (25) 


^ [9a^ 7 ^ 3 ( 0 ) -|-4a^Ar2(a) +a^Ki{a) 

+2b {ia'^K2{a) + a^Ki{a) - ia^K^i^a) + a^K2{a)) 
6 ^ {3a^K2{a) + a^Ki{a))] (26) 
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Fig. 4. The ratio of the particle density calculated to first 
order in {q — 1) normalized to the particle density of a Tsallis 
gas as a function of the temperature for different values of the 
parameter q. The mass is taken as being the pion mass. The 
values of the parameter q are 1.01 for the dashed line, 1.08 for 
the dot-(long)dashed line, 1.1 for the dot-(short)dashed line 
and 1.15 for the dotted line. 



0.01 but rapidly deviates from the full Tsallis distribu¬ 
tion for larger values of q. Also here, for {q — 1) « 0.1 
the deviations are of the order of 20 % as can be seen 
from Fig. For comparison we show in Fig. the first 



Fig. 6. The ratio of the energy density calculated to first 
order in (q — 1) normalized to the energy density of a Tsallis 
gas as a function of the temperature for different values of the 
parameter q. The mass is taken as being the pion mass. The 
values of the parameter q are 1.01 for the dashed line, 1.08 for 
the dot-(long)dashed line, 1.1 for the dot-(short)dashed line 
and 1.15 for the dotted line. 


order expansion compared to the Boltzmann expression, 
(e^ + (9 — l)n^)/e'®, as a function of the temperature T 
for several values of the parameter q. As in the previous 
case the deviations are most pronounced for small values 
of the temperature. Finally, the pressure is given by 

+ {q- 1)P^ (27) 


pB 


ge^T^a^ 7^2(0) 
2tt'^ 


(28) 


Fig. 5. The ratio of the particle density calculated to first 
order in {q — 1) normalized to the particle density of a Boltz¬ 
mann gas as a function of the temperature for different values 
of the parameter q. The mass is taken as being the pion mass. 
The values of the parameter q are 1.01 for the dashed line, 1.08 
for the dot-(long)dashed line, 1.1 for the dot-(short)dashed line 
and 1.15 for the dotted line. 


In Fig. we show the ratio of the energy density to 
first order in (9 — 1 ) to the full energy density as given 
by the Tsallis distribution, (e^ + (9 — l)e^)/e for several 
values of 9 indicated in the figure as a function of the tem¬ 
perature T. 

Again, as noted previously for the particle density, it can 
be seen that the expansion in ( 9 — 1 ) is excellent if ( 9 — 1 ) = 


P^ = [a^K^ia) + 803 ^ 3 ( 0 ) 

-2a%K3{a) + a%'^K2[a) + 2a^hK2{a)\ (29) 

In Fig. [^we show the ratio of the pressure to first or¬ 
der in (9 — I) to the full pressure as given by the Tsallis 
distribution, {P^ + (9 — l)P^)/P for several values of 9 
indicated in the figure as a function of the temperature T. 
Again, as noted previously, it can be seen that the expan¬ 
sion in (9 — 1 ) is excellent if (9 — 1 ) = 0.01 but rapidly 
deviates from the full Tsallis distribution for larger values 
of 9 . Also here, for (9 — 1 ) « 0.1 the deviations are of the 
order of 20 % as can be seen from Fig. Again, we show 
in Fig. the first order expansion compared to the Boltz¬ 
mann expression, (P^ + (9 — ^)P^)/P^: as a function of 
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Fig. 7. The ratio of the energy density calculated to first order 
in (q— 1) normalized to the energy density of a Boltzmann gas 
as a function of the temperature for different values of the 
parameter q. 



Fig. 8. The ratio of the pressure calculated to first order in 
(g — 1) normalized to the pressure of a Tsallis gas as a function 
of the temperature for different values of the parameter q. 


the temperature T for several values of the parameter q. 
In this case the deviations are most pronounced for small 
values of the temperature. 

The entropy density can be obtained from the above 
expressions by using the thermodynamic relation 

€ + P = Ts + fj.n (30) 

and is not being shown here. 

Before closing the section, we comment about the va¬ 
lidity of the present expansion up to the second order in 
(g — 1). As we have already seen, for truncation of the ex¬ 
pansion up to the first order with pure Tsallis distribution. 



Fig. 9. The ratio of the pressure calculated to first order in 
(g — 1) normalized to the pressure of a Boltzmann gas as a func¬ 
tion of the temperature for different values of the parameter g. 


one needs to satisfy two conditions, i.e. 


|l-9lf <1 


(31) 


and 


|l-<zl(fl <2 


(32) 


With the thermodynamically consistent Tsallis distri¬ 
bution, the second condition turns out to be \q{l—q)\{E/T)'^ 
< 2. For expansion up to second order in (g — 1) with 
modified Tsallis distribution, the condition becomes g^|l — 
g|(A/T)3 <3. 

Given the two values (the highest and the lowest) of 
(g — 1) used in the present analysis, we want to put an 
upper bound in E jT until which the expansion will be re¬ 
liable. If (g — 1) = 0.15, Eq.(31) gives E/T < 6.7 and 
the modified condition for expansion up to first order 
gives {E/T)"^ < 11.59. That means for T = 0.1 GeV, 


the expansion will be reliable up to E 
(g- 1) = 0.15. 


0.3 GeV when 


If (g - 1) = 0.01, Eq.@ gives E/T < 100 and the 
modified condition for expansion up to first order gives 
(E/T)'^ < 198. That means for T = 0.1 GeV, the expan¬ 
sion will be reliable up to A « 1.4 GeV when (g — 1) = 
0.01. The permissible values of (g — 1) and E/T for reli¬ 
able expansion up to the second order in (g — 1) are shown 
to reside inside the smaller area filled with slanted lines. 
This is shown in Fig. [TOl 

Hence, we conclude from the above discussion that the 
smaller the g value the more reliable the expansion be¬ 
comes. 
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Fig. 10. (colour online) The region of validity for the expan¬ 
sion in (g — 1) up to second order. The area under the solid line 
denotes the region where only the condition \1 — q\E/T < 1 is 
satisfied. The common overlapped area (i.e. the area under the 
dotted line) depicts the region within which all the conditions 
for expansion up to second order are satisfied. 


5 Inclusion of Flow to First Order in (g — 1) 


In order to see how the inclusion of flow could improve 
the description of the transverse momentum distributions 
obtained in Pb-Pb collisions, we have included a con¬ 
stant flow velocity, /3. Assuming space-like freeze-out sur¬ 
face, the invariant yield is given by (see Appendix for the 
derivation). 

1 dN _ gV 
Pt dprdy (27r)2 

^‘2T[rIo{s)Ki{r) - sIi{s)Ko{r)] 

-{q - l)Tr‘^Io{s)[Ko{r) + K 2 ir)] 

+4:{q— 1) TrsIi{s)Ki{r) 
-iq-l)Ts^Ko{r)[Io{s) + his)] 

+ ^^^rr3jo(s)[if3(r) + 3i^i(r)] 
_ ^(^-^) Trh[K 2 ir) + i^o(r)]/i(s) 

+ ^(^~l) rs^r[/o(s) + /2(s)]i^i(r) 

- ^^^Ts3[/3(s) + 3/i(s)]XoW| (33) 

where 


^rriT 

T 

^VPT 

T 


(34) 

(35) 


In{s) and Kn{r) are the modified Bessel functions of the 
first and second kind. Now, in this formula, the freeze- 
out surface has been considered to be space-like and so 
the integration over the freeze-out surface turns out to 
be trivial. For a more detailed treatment of the freeze-out 
surface in this context, the readers are referred to Ref. 

El- 



Fig. 11. Fits to the normalized differential tt yields as mea¬ 
sured by the ALICE collaboration in (0 — 5)% Pb-Pb collisions 
at y'SNN = 2.76 TeV [10]. The fit with the Tsallis distribution 
including flow keeping terms to first order in [q — 1) (dashed 
line). The flow velocity is fixed at /3 = 0.609, with T = 146 
MeV, q = 1.030 and the radius of the volume is R — 29.8 
fm. The solid line is the Tsallis distribution without flow as 
given in Fig. The lower part of the figure shows the differ¬ 
ence between model (M), i.e. Tsallis with flow up to first order 
in {q — 1), and experiment (E) normalized to the model (M) 
values. 


The comparison between model and experiment is quite 
good with notable deviations at small values of the trans¬ 
verse momentum pT and again above values of 2.5 GeV 
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(see Fig. 11). These could easily be attributed to the coarse This is an elementary integral over a polynomial func- 
way of treating transverse flow. More detailed investiga- tion. The result is 
tions have been carried out in 1221. 


6 Summary and Conclusion 

The Tsallis distribution describes extremely well the trans¬ 
verse momentum distributions in p - p collisions at high 
energies. All fits performed so far show that the parameter 
q is always close to 1. In view of this, we have presented 
in this paper a series expansion of quantities relevant in 
the analysis of high energy physics in (g — 1). The Tsal¬ 
lis distribution itself has been obtained to second order 
in (g — 1). A rough comparison with experimental data 
has been done for the transverse momentum distributions 
obtained in p-p collisions. In the case of Pb-Pb collisions 
we have given an estimate using flow with a fixed flow ve¬ 
locity. In most cases the series expansion turns out not to 
be a useful description of transverse momentum distribu¬ 
tions but it could be useful in analyses where a comparison 
and detailed investigation is needed when comparing the 
Tsallis and the Boltzmann distributions. Furthermore, a 
systematic study of the identified particle yield due to p—p 
and heavy-ion collision at RHIC and LHC has been done 
in [35] where the flow formula (Eq. (33)) has been used to 
fit the spectra. 


dN 

dy 


y=o 


gv 

(2^)2 


l + (9-l) 

2^3 


m — fi 
T 


-2 


(2g-3)(g-2) 

- 


2-(g-2) 


m — p 
T 


+ 2 


m — p 


T 


g-2 


1 + 


m — p 
T 


+ Tp^ 


This makes it possible to express the volume V in terms 


dy 


y=o 


from Eq. (13) in the text. Replacing the volume 


term in Eq. (12) in terms of Eq. (39), we get Eq. (14). 


B Expansion in (g — 1) of Tsallis distribution 


(39) 
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Appendix 


E - p 

q-l = x\ = <P 

l + {q- = l+x<P = f{x) 


q 

■g-I 


I -I- a; 


= g{x) 


(40) 

(41) 

(42) 


A Derivation of Eq. (13) and Eq. (14) 


The result quoted in Eq. ( [T^ can be obtained by a change 
of variable from px to x which is defined by 


The relevant distribution function can then be written as 


= l + (g-l) 


rriT — ^ 


T 


(36) 


This leads to the following integral 

dN 
dy 

gV 


y=0 


with 


/■°° , dN 
lo dprdy 


(2^)2 

gV 


y=o 

dpT Ptitt-t 


1 + (<Z-1) 


rriT — ^ 


T 


(27r)2 g - 1 



X — 1 

/ dx X 9-1 

J a 

- tT + P 

Li-1 J 


a = 1 -I- (g — 1) 


m — II 


l + (9-l) 


E — p 
T 


'9-1 l+x 

= [1 -I- 

= /(x)®(") 

= E{x) 


(43) 


Expanding E{x) in Taylor series about a; = 0, we get 


E{x = 0) = e 


<2> 


(37) 

(38) 


dE{x) 


dx 


x—0 


= -{-2 + d>)cPe-^ 


(44) 

(45) 


^ (24 - 20<? + 3<?>2) e-^ (46) 

dx^ x^o 12 ^ ' 
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for 0{x^), 0(x) and 0{x‘^) respectively. The final results 
can be summarized in the following equation 


And hence, 

B , g{q-l)eTT^ 


n = n -I- 


l + (9-l) 


E — g, 
T 


dTT^ 


du! u) {uJ^ — a^) ^ 


E-fJ 

~ e T 


/ ^ E — g\ e-m 

+ ( 9 -I- 2 +I e ^ 


2 T 


T 


+ 


{q-i? 1 
2 ! 12 


E — g 
T 


E — u, f E — u 

24 - 20- - + 3 - - 

T \ T 


e T 


(47) 

(48) 


For completeness we also quote the following result 

E — g 


(— 2 a; 2h + uj'^ -\-— 2ujb) e 
n = n^ + g(g-l)e»r3 _ 2a^K,{a) 

—40^6X2(0) + 'ia^K^{a) + a^K 2 {a) + a'^b'^K 2 {a) 

-2a^bKi{a)] (52) 


Here the following form of the modified Bessel function of 
second kind |45| is used: 


on —1/„ -|\| roo 

= ( 2 n- 2 )!a» A ^ (^3) 

On the other hand, by using Boltzmann Statistics one 
obtains. 


l + (g-l)- 


T 


_ E — fj. 

~ e 


n^ = 


ge^T^a^K2{a) 

2'k'^ 


(54) 


2 V T 


+ ( 9 - 1)0 ( I e E 


E — g 


{q-ir 1 

E — g. 

2 ! 12 

T 

. f E — u\ 

T 

-j +3 


C.2 Pressure 

To first order in [q— 1) the pressure is given by: 


T 


E — fj. 


P = 


j3 I" 

(2^)3 7 ^ ^ SA 


9 


e ^+( 5 - 1 ) 


E — g 
2T 


(49) 


„ E — n\ E-u 

- 2 +^^le ^ 


(55) 


C Tsallis Thermodynamics 

C.l Particle number density n 

Up to 0{q — 1), the number density can be written as: 

E — fJL 


Using the same definitions as in Eq. (511, we get 

, g{q-l)e^T^ 


P = + 


127r2 


J doj {uj"^ — a^) ^ 


(—2a; -I- 26 -I- a;^ -I- 5^ — 2ujb) e 
^pB ^ g{q-l)e^T^ [04^2(0) + ia^K^ia) 


-'‘+(9-1)- 
-2+^\ e-^ 


2T 


T 


(50) 


dTT^ 

—2a^bK3{a) + a^6^ 7 ^ 2 ( 0 ) + 2a^bK2{a)] (56) 

Using Eq. (51) along with Eq. ( [5^ as well as another 
representation of the modified Bessel function, 

on„| roo 

= (57) 


We define. 


Eq. (53) can be obtained from Eq. (57) by dint of partial 


integration. The Boltzmann pressure density is, similarly 
given by: 
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C.3 Energy density e 


D Momentum Distribution 


e = 


(27r)3 


(fp E 


e-m , — n 

e E +{q-l)- 


2T 


-2 + 


E — fi\ _ JEnii 


T 


e E 


(59) 


Using Eq. (511 along with Eq. (53), 

e = -I- ^ [ 9 a^-ftr 3 (a) -I- 4 a"‘itr 2 (a) + a^Ki{a) 

+2b ( 80 ^ 1 ^ 2 ( 0 ) -I- a^Ki{a) — Za^K^{a) + a^K2{a)) 


(ia'^K2{a) -\- a^Kl{a))^^ 


with 


3 " {3a^K2{a) + a^Ki{a)) 


2 Tr^ 


C.4 Checking Thermodynamic Consistency 


The invariant particle yield is given by, 

E — p 


dN 


l + (9-l)- 


T 


(67) 


Where C = . Assuming g — 1 «C 1 we are to expand 

it in Taylor series. Let, 


dN 


pxdpTdydcf) 


= CE 




+ ^^(-2 + <?>) 




(60) 




+ ^j^{2A-20d> + 3cP^)e-^ + 


( 68 ) 


(■ 02 ) where y is rapidity and (j) is the azimuthal angle of emis¬ 
sion. Now, parameterizing energy E in terms of rapidity 
y i.e. putting E = ruTCOshy and integrating over y we get 

Eq. (figl). 


We know that for thermodynamic variables n and P at 
non-zero /i we have. 


E {q- l)n^ -b 


n = dP/dpL 
.. = dP^Idp P 

(q- l)5PV5/i- 


(62) 


■ (63) 

where the superscripts denote the order of {q — 1 ) in the 
expansion. Since the terms in the expa nsio n are linearly 
independent, the identity given by Eq. (62) is to be sat¬ 


isfied at every order of {q— 1). Since our expansion is up 
to C>(g — 1 ), we will check the identity for that order now. 
Rearranging and reordering the 0{q—l) terms in Eq. (52) 
we get. 


[a^A: 3 (a) -b a'‘A: 2 (a) 

—Aa^bK2{a) + 2a^b{K^{a) — Ali(a)) 
—2a^bK^{a) -b a^b^K 2 {a) — 6 a^K 2 {a) 

+ 2a^iK3{a) - Kiia))] 

(64) 

Using the recursion relation for the modified Bessel’s func¬ 
tions, 


‘2i71 

Kn+i{a) - Kn-i{a) = —Ar„(a) 


(65) 


[a^K^ia) + a^K 2 {a) 

+Aa^bK2{a) — 2a?bK3[a) + a^h^K2{a) 

+ 2 a'^K 2 {a)\ 

dP^ 


dy 

Hence, proved. 


( 66 ) 


E Flow in Tsallis Distribution 

We use the following ansatz (in cylindrical polar coordi¬ 
nates) for introducing flow inside our calculations: 


= {mTCOshy,pTCOs4>,pTsi'a(j),mTSuihy) (69) 
= (qcosh^, yucosa, yusina, ysinhC) (70) 

where {Qy is the (space-time)rapidity of particles (fluid- 
element) and V is the velocity of fluid. Now, to include 
flow inside the Tsallis distribution, we replace E p^u^ 
assuming temperature to be scalar. The dot product then 
becomes. 


p^Ufj, = ymTCOsh(y — C) ~ ^vpTCOs{(j) — a) 


(71) 


Now putting Eq. ( |7T| ) in Eq. (18) up to 0{q — 1) and 
integrating over </> and C, we get. 


I dN _ gV 
Pt dprdy (2-7r)2 

{2T[r/o(s)A:i(r) - s/i(s)A:o(t’)] 

-{q - l)Tr'^Io{s)[Ko{r) + K 2 {r)] 

-b4(q — 1) TrsIi{s)Ki{r) 
-{q-l)Ts^K3{r)[Io{s)pl2{s)] 

+ ^l^Tr^I^{s)[K 3 {r) + 3ifi(r)] 
-^l^Tr^^[K2{r) + ifo(0]A(s) 

+ /2(s)]77i(r) 

- ii^rs3[/3(s) + 3Ji(s)]7^o(r)| (72) 
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where 


'yrriT 

r = - 

T 

JVPT 

® ^ T 

In{s) is the modified Bessel function of first kind. 
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